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Abstract 

We consider a fading AWGN 2-user 2-hop network where the channel coefficients are independent and identically 
\ distributed (i.i.d.) drawn from a continuous distribution and vary over time. For a broad class of channel distributions, 

we characterize the ergodic sum capacity to within a constant number of bits/sec/Hz, independent of signal-to-noise 
ratio. The achievability follows from the analysis of an interference neutralization scheme where the relays are 
partitioned into M pairs, and interference is neutralized separately by each pair of relays. When M — 1, the 
proposed ergodic interference neutralization characterizes the ergodic sum capacity to within 4 bits/sec/Hz for i.i.d. 
uniform phase fading and approximately 4.7 bits/sec/Hz for i.i.d. Rayleigh fading. We further show that this gap 
can be tightened to 41og7r — 4 bits/sec/Hz (approximately 2.6) for i.i.d. uniform phase fading and 4 — 41og(^) 
bits/sec/Hz (approximately 3.1) for i.i.d. Rayleigh fading in the limit of large M0 

Index Terms 

Amplify-and-forward, approximate capacity, ergodic capacity, fading, interference neutralization, two unicast, 
Q . two-user two-hop networks. 

I. Introduction 

>. 

CN . In recent years, there has been significant progress towards understanding fundamentals of multi-source single- 
hop networks CD-El H Following up on these successes for single -hop networks, more recent and emerging work 
■ has considered multi-source multi-hop networks El-Hi. For multi-source multi-hop networks, interference can be 
cancelled by aligning multiple paths through the network, a technique referred to as interference neutralization. 
Proper exploitation of such interference neutralization is the key for an approximate capacity and the optimal 
<N i de grees of freedom (DoF) characterization [|3j — QSj . Recently, for 2-user 2-hop networks, interference neutralization 
. . . combining with symbol extension was used to show that two relays suffice to achieve the optimal DoF J5j. In spite 
. £h ! °f recent progress in this area, the best known capacity characterization for fully connected 2-user 2-hop networks 
^ ■ is to within o(log(SNR)) bits/sec/Hz |5 ], which can be arbitrarily large as the signal-to-noise ratio (SNR) increases. 

The aim of this paper is to tighten the capacity gap of 2-user 2-hop networks to within a constant number 
of bits/sec/Hz, independent of SNR. Our achievability is based on ergodic interference neutralization (6), which 
is similar to ergodic interference alignment [3] applied to multi-source single-hop networks. Suppose that the 
sources transmit their signals at time t through the first-hop channel matrix H[i]. Then the relays amplify and 
forward their received signals with an appropriate delay r through the second-hop channel matrix G[t + t] such 
that G[i + r]H[i] becomes an approximately diagonal matrix with non-zero diagonal elements. This approach can 
completely neutralize interference in the finite SNR regime. 

Assuming independent and identically distributed (i.i.d.) channel coefficients, the proposed ergodic interference 
neutralization characterizes the ergodic sum capacity to within a constant number of bits/sec/Hz for a broad class 
of channel distributions. For instance, when the number of relays L is equal to two, it achieves the ergodic sum 
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TABLE I 

New approximate capacity results and the existing DoF and approximate capacity results. 





A"-user 
interference channel 


2-user 2-hop network 
with 2 relays 


A"-user 2-hop network 
with L relays 


A" -user if -hop network 
with AT relays at each layer 


DoF 


f m 


2 (5) 


Generally unknown 
K if L > K(K - 1) + 1 (9) 
L if K -> oo (101 


Generally unknown 
K for isotropic fading \6\ 


Ergodic capacity 
for uniform 
phase fading 


Exact capacity 1111 


4 bits/sec/Hz gap 


Generally unknown 
2.6 bits/sec/Hz gap 
if K = 2 and L — > oo 


Unknown 


Ergodic capacity 
for Rayleigh fading 


1.3K bits/sec/Hz gap 


4.7 bits/sec/Hz gap 


Generally unknown 
3.1 bits/sec/Hz gap 
if K = 2 and L — > oo 


Unknown 



capacity to within 4 bits/sec/Hz for uniform phase fading and approximately 4.7 bits/sec/Hz for Rayleigh fading. 
As L increases, we narrow the corresponding gap in our analysis. Specifically, this gap is given as 4 log 7r — 4 
bits/sec/Hz (approximately 2.6) for i.i.d. uniform phase fading and 4 — 41og(^) bits/sec/Hz (approximately 3.1) 
for i.i.d. Rayleigh fading in the limit of large L. We also notice that a similar analysis is applicable for the 
if -user interference channel and show that ergodic interference alignment in (3] characterizes the ergodic sum 
capacity assuming that all sources employ uniform power allocation across time to within (| log 6) if bits/sec/Hz 
(approximately 1.3if) for i.i.d. Rayleigh fading. Table U summarizes the new approximate ergodic capacity results 
of this paper and the existing DoF and approximate capacity results. 

A. Related Work 

1 ) Degrees of freedom: In seminal work |2 ], interference alignment has been proposed to achieve the optimal DoF 
of the if -user interference channel with time-varying channel coefficients. The concept of this signal space alignment 
has been successfully adapted to various network environments, e.g., see ifTOl . lfl2l - lfl"8l and the references therein. 
It was shown in |fl~9l , 1201 that interference alignment can also be attained on fixed (not time-varying) interference 
channels. 

In spite of recent achievements on interference channels or multi-source single-hop networks, understanding of 
multi-source multi-hop networks is still in progress. The work |5 ] has exploited interference alignment to neutralize 
interference at final destinations, which is referred to as aligned interference neutralization, and showed that the 
optimal 2 DoF is achievable for 2-user 2-hop networks with 2 relays. This result has been recently generalized 
to two unicast networks 0, (H. For more than two unicast, the optimal DoF is in general unknown except for a 
certain class of networks. For the if -user 2-hop network with L relays, interference can be completely neutralized 
if L > if (if — 1) + 1 [9]. Similar concept of ergodic interference alignment has been proposed for interference 
neutralization in showing that ergodic interference neutralization achieves the optimal DoF of if -user if -hop 
isotropic fading networks with K relays in each layer. 

2 ) Beyond degrees of freedom: The DoF discussed previously is a fundamental metric of multi-source networks 
especially for high SNR, which characterizes capacity to within o(log SNR) bits/sec/Hz. Depending on the opera- 
tional regime, however, the gap of o(log SNR) bits/sec/Hz in practice can be significant and achieving the optimal 
DoF may not be enough. For the 2-user interference channel, for instance, time-sharing between the two users 
can also achieve the optimal one DoF. On the other hand, a simple Han-Kobayashi scheme can tighten the gap 
to within one bit/sec/Hz 100, which provides an arbitrarily larger rate compared with the time-sharing for a certain 
operational regime and channel parameters. Consequently, several works have recently established tighter bounds 
on the gap from capacity (H, 11211 - 11261 to provide a universal performance guarantee, independent of SNR and 
channel parameters. 

A similar flavor of such bounds on the gap from capacity concerns time-varying channel models. The recently 
proposed ergodic interference alignment in (3J makes interference aligned in the finite SNR regime and, as a result, 
provides significant rate improvement compared with the conventional time-sharing strategy in the finite SNR 
regime. Ergodic interference alignment was shown to achieve the ergodic sum capacity of the if -user interference 
channel for i.i.d. uniform phase fading 0. For the if -user finite field interference channel (with time- varying 
channel coefficients), the idea of ergodic interference alignment was independently proposed by Nazer et al. ifTTTl 
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TABLE II 

Summary of notation 



A J ( or a J ) 


Transpose of A ( or a) 


A f ( or a T ) 


Conjugate transpose of A ( or a) 


det(A) 


Determinant of A 


I 


Identity matrix 


J 


V-i 


re(a)( or im(a)) 


Real (or imaginary) part of a 


\a\ 


Absolute value of a 


a* 


Complex conjugate of a 


[a] 


Floor of a ([a\ — max{i g Z\x < a}) 


card (.4) 


Cardinality of A 




Circularly symmetric complex Gaussian distribution with mean /i and variance o* 



and Jeon and Chung E71 in two slightly different versions. In |[28l . ergodic channel pairing was applied to tighten 
the gap from the ergodic capacity for fading multihop networks showing a gap depending only on the number of 
nodes in a layer, instead of the total number of nodes in a network. 



B. Paper Organization 

The rest of the paper is organized as follows. In Section |II1 we introduce the fading 2-user 2-hop network model 
considered in this paper and formally define its ergodic sum capacity. In Section [Till we first state the main results of 
the paper, approximate ergodic sum capacities of fading 2-user 2-hop networks. In Section [IVJ we explain ergodic 
interference neutralization and its achievable rate. In Section |VJ we prove the approximate ergodic sum capacity 
results in Section [III] based on the achievability in Section |IV] Finally, we conclude in Section [Vl] and refer some 
technical proofs to the appendices. 

II. Problem Formulation 

In this section, we explain our network model and define its sum capacity. Throughout the paper, we will use A, 
a, and A to denote a matrix, vector, and set, respectively. The notation used in the paper is summarized in Table 

in 



A. Fading 2-User 2-Hop Networks 

We study the 2-user 2-hop network depicted in Fig. Q] in which each source wishes to transmit an independent 
message to its destination with the help of L relays, where L > 2. The input-output relation of the first hop at 
time t is given by 

y R [t]=H[t]x[t] + z R [i\, (1) 



where 



H[t] 



h 2 ,i[t) 



hi, 2 [t] 

h2,2 [t] 



(2) 



h L M h L , 2 [t] 

is the L x 2 dimensional complex channel matrix of the first hop at time t, y_R[i] = [y_R,i[t], ■ ■ ■ ,yR,h[t]] T is the 
L x 1 dimensional received signal vector of the relays at time t, x[t] = [xi[t], X2[t]] T is the 2x1 dimensional 
transmit signal vector of the sources at time t, and z#[i] = [^H,i[i]> • • • j "2r,l[£]] T is the L x 1 dimensional noise 
vector of the relays at time t. Similarly, the input-output relation of the second hop at time t is given by 



y[t} = G[t]x R [t]+z[t], 



where 



G[t] 



9i,i [t] 91,2 [t] 

92,1 [t] 92,2[t] 



9l,L[t] 
92 At] 



(3) 



(4) 



4 



Wi-* Si 



Xi[t] 



x 2 [t] 



H[t] 



Z R,1 M 



vrAA 



Hi 



XR,l[t] 



XR,2 [t] 



XR,h[i\ 



G[t] 



zrAA 

Fig. 1. Gaussian 2-user 2-hop network with L relays. 



Zl [t] 
1 



2/2 [t] 
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is the the 2 x L dimensional complex channel matrix of the second hop at time t, y[t] = [yi[t],y2[t]] T is the 
2x1 dimensional received signal vector of the destinations at time t, x#[t] = [xr,i[£], • • • , XR^ft]] 7 is the L x 1 
dimensional transmit signal vector of the relays at time t, and z[t] = [z\[t], Z2[t]] T is the 2x1 dimensional noise 
vector of the destinations at time t. We assume that the elements of z#[i] and z[t] are i.i.d. drawn from £A/"(0, 1). 
Each source and relay should satisfy the average power constraint P, i.e., E[|xj[t]| 2 ] < P for i G {1,2} and 
E[\x Rd [t]\ 2 ]<P for j e {!,■■■ ,L}. 

We assume that channel coefficients are i.i.d. drawn from a continuous distribution f(x), x 
independently over time. Without loss of generality, we assume that E[|/ijj [t]| 2 ] = 1 and Eflg^j] 
i G {l,--- ,L} and j G {1,2}. We further assume that the sources do not know any channel state information 
(CSI) and the relays and the destinations know global CSI. That is, at time t, each relay and destination knows 
H[t] and G[t}. 



G C, and vary 
]\ 2 } = 1 for all 



B. Ergodic Sum Capacity 

Based on the network model, we consider a set of length-n block codes. Let W% be the message of source 
i uniformly distributed over {1, • • • ,2 nR '}, where R{ is the rate of source i. A rate pair (R\,R2) is said to be 
achievable if there exists a sequence of (2 nRl , 2 nR2 , n) codes such that the probabilities of error for W\ and W% 
converge to zero as n increases. Notice that since channel coefficients are i.i.d. varying over time, an achievable 
rate pair (Jf2i, JR2) is given as in the ergodic sense, i.e., the expectation over random channel coefficients. The 
ergodic sum capacity C sum is defined as the maximum achievable ergodic sum rate. Unless otherwise specified, 
an achievable sum rate or the sum capacity in this paper mean an achievable ergodic sum rate or the ergodic sum 
capacity, respectively. 



III. Main Results 

In this section, we first introduce our main results. Let M := |_§J- As will be explained in Section |IVJ we only 
use 2M relays among the total number L of relays for the achievability. That is, the achievability is based on an 
even number of relays. Without loss of generality, we assume that relay 1 to relay 2M are used for relaying. The 
achievability follows from ergodic interference neutralization based on amplify-and-forward relaying in which 2M 
relays are partitioned into M pairs and interference is neutralized separately by each pair of relays. In order to 
describe the proposed ergodic interference neutralization and its achievable sum rate, for m G {1, • • • ,M}, we 
denote 

H m [t] 

and 

G m [t] 



^2m-l,l[i] 
^2m,l [t] 



h2m-l,2\t\ 
h2m,2 [t] 



9l,2m-l[t] <?2,2m-lM 
9l,2m[t] 92,2m\P] 



(5) 



(6) 
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which are the 2x2 dimensional channel matrices at time t from the sources to relays 2m 
relays 2m — 1 and 2m to the destinations, respectively. 



1 and 2m and from 



A. Achievable Sum Rate 

The following theorem states an achievable symmetric rate of the fading 2-user 2-hop network. 

Theorem 1: For the fading 2-user 2-hop network with L relays, 



log 



1 + 



Pi 2 E"=i|det(H r 



V 



l + o- 



AFA 



(V) 



P n 2 
1+2P' AF, 



7 2 E^=i(l^( 



m— 1)+3— i,3— i I 



l + \h 



2(m-l)+i,3-i| J> 



is achievable for i £ {1,2}, where M = L§J » 7 
and the expectation is over the channel coefficients. 

Proof: The proof is in Section [IVJ ■ 

The most important aspect is that there is no residual interference after ergodic interference neutralization, 
meaning that interference can completely be neutralized at finite SNR. Moreover, from the block-wise coherent 
combining gain shown as (Em=i I det(H m )|) 2 in (O, the received signal power increases as the number of pairs 
M increases. Although there is noise amplification due to amplify-and-forward relaying given as a\-p i in d7), this 
additional noise results in a constant number of bits/sec/Hz loss for a broad class of channel distributions, which 
will be proved in Section [V] 

For notational convenience, let 

/ w 
log 1 + 



i=l 



^7 2 (E*=i|det(H, 



V 



1 + a 2 



AFA 



(8) 



which is the achievable sum rate from Theorem [TJ For comparison, we consider the ergodic capacity of the multiple - 
input multiple-output (MIMO) channel from the sources to the relays, that is 



Rn 



logdet(I + PHH j ; 



(9) 
fl29l 



Since the channel coefficients are i.i.d. and the sources do not know CSI, C sum is upper bounded by R m 
The following example illustrates R m and -R m i mo for i.i.d. Rayleigh fading, i.e., f(x) follows CA/"(0, 1). 

Example 1 (Sum rate: Rayleigh fading): Figure |2] plots R m and R mimo for i.i.d. Rayleigh fading. Two important 
aspects can be observed in the figure. First, for a fixed number of relays L, the sum rate gap i? m i mo — Rm appears to 
be upper bounded by some constant independent of power P, which suggests that the proposed ergodic interference 
neutralization can achieve the ergodic sum capacity to within a constant number of bits/sec/Hz independent of P. 
Second, for a fixed P, the sum rate gap i? m i mo — R m appears to decrease with increasing L, which suggests that 
this approximate capacity characterization can be tightened as the number of relays L increases. 

Both observations in Example Q] are established in this paper and shown to hold beyond the case of Rayleigh 
fading for any fading model for which f(x) is only a function of The following two subsections describe our 
approximate capacity results characterizing the ergodic sum capacity to within a constant number of bits/sec/Hz, 
independent of P. 



B. Approximate Ergodic Sum Capacity for L = 2 

In this subsection, we assume L = 2. We first consider i.i.d. uniform phase fading in which hij[t] = exp(j#j j[t]) 
and gj t i[t] = exp(jipj t i[i\), where 0ij[t] and <fj,i\t\ are uniformly distributed over [0, 2tt) for all i,j G {1,2}. 
Although uniform phase fading violates the channel assumption in Section III- A I i.e., f(x) is continuous over 
x G C, we can slightly modify the proposed ergodic interference neutralization and show that Theorem Q] still 
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Fig. 2. The achievable sum rate Ri n and its upper bound i? m i mo for i.i.d. Rayleigh fading when L — 2,4, 16, 64. 



holds. The detailed modification is given in Appendix I. The following theorem characterizes an approximate 
ergodic sum capacity for i.i.d. uniform phase fading. 



Theorem 2: Consider the fading 2-user 2-hop network with L = 2 relays. If hij[t] = exp(j#jj[i]) and gji[t] 
exp(j(pj t i[t]), where Oij[t] and <pjj[t] are uniformly distributed over [0, 2ir) for all i,j G {1,2}, then 



Rin <4 



(10) 



for any P > 0. 

Proof: The proof is in Section IV-AI ■ 

Example 2 (Gap for L = 2: Uniform phase fading): Figure [3] plots i? m imo — ^in with respect to P for i.i.d. 
uniform phase fading (the closed forms of i? m i mo and R[ n are given by (T37T ) and (1381 ). respectively). As proved by 
Theorem [2l the proposed ergodic interference neutralization achieves C sum to within 4 bits/sec/Hz for i.i.d. uniform 
phase fading. This theoretical gap coincides with the actual gap -R m i mo — Ri n at high SNR, i.e., limp^ > . 00 {R m i mo — 
Rin} =4. 

Based on the bounding techniques used in proving Theorem |2j we characterize an approximate ergodic sum 
capacity for a class of channel distributions satisfying that f(x) is only a function of \x\. Specifically, for a given 
set of channel amplitudes, we first upper bound the gap i? m i mo — -Rin by averaging out the effect of phase fading. 
Then we further apply additional bounding techniques to obtain an upper bound, independent of power P. 

Theorem 3: Consider the fading 2-user 2-hop network with L = 2 relays. If f(x) is only a function of then 



R\„ < 2 E 



log 



B(A + yjA 2 - G 2 ) 



+ 2 



(11) 



for any P > 0, where 



A = |/ii,i| 2 |/i 2 , 2 | 2 + \hi, 2 \ 2 \h 2 ,i\ 2 , 
B = \h hl \ 2 + \h 2 ,i\ 2 + 2, 

G = 2|/ii 1 1 1/112I |^2,i 1 1^<2,2|) 



(12) 
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Fig. 3. Gap from the sum capacity for i.i.d. uniform phase fading when L — 2. 
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Fig. 4. Gap from the sum capacity for i.i.d. Rayleigh fading when L = 2. 



and the expectation is over the channel coefficients. 

Proof: The proof is in Section IV-Ai ■ 

The presented gap in Theorem [3] only depends on the amplitude distribution of channel coefficients, which 
provides universal performance guarantee regardless of power P. The following example evaluates the presented 
gap for i.i.d. Rayleigh fading. 

Example 3 (Gap for L = 2: Rayleigh fading): Figure @] plots -R m i mo — R m with respect to P and also plots its 
upper bound in Theorem [3] for i.i.d. Rayleigh fading. Since there is no closed form, we evaluate the bound in 
Theorem [3] by simulation, which approximately provides 4.7 bits/sec/Hz gap. Simulation result shows that the 
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b percent of the ergodic sum capacity at SNR 



C. Approximate Ergodic Sum Capacity as L — > oo 

In this subsection, we focus on an approximate ergodic sum capacity as the number L of relays increases. Again, 
we first consider i.i.d. uniform phase fading and then consider a class of channel distributions satisfying that f(x) 
is only a function of \x\. 

Theorem 4: Consider the fading 2-user 2-hop network with L relays. If hij[t] = exp{j9ij[t]), gji[t] = &xp{j<Pj,i[t]), 
and 0ij[t] and tpj,i[t] are uniformly distributed over [0, 2tt) for all i G {1, • • ■ , L} and j G {1, 2}, then 



lim {C sum 

L— >oo 



Rin} <41og7T-4 



(13) 



for any P > 0. 

Proof: The proof is in Section IV-BI ■ 

Example 4 (Gap as L — >• oo: Uniform phase fading): Figure [5] plots the gap -R m i mo — -Rin for i.i.d. uniform phase 
fading with respect to L. As shown in the figure, this gap decreases as L increases and eventually converges to 
41og7r — 4 (approximately 2.6) regardless of P, which was proved in Theorem 01 Therefore the proposed ergodic 
interference neutralization characterizes C sum to within 4 log n — 4 bits/sec/Hz in the limit of large L. Compared 
to 4 bits/sec/Hz, the sum capacity gap for L = 2 in Theorem the result shows that the sum capacity gap can be 
tightened as L increases. 



Theorem 5: Consider the fading 2-user 2-hop network with L relays. If f(x) is only a function of then 

lim {C sum - R in } < 4 - 41og (E[| det(H x )|]) (14) 

for any P > 0. 

Proof: The proof is in Section IV-BI ■ 

Example 5 (Gap as L — >• oo: Rayleigh fading): Figure [5] plots i? m i mo — -Rin for i.i.d. Rayleigh fading with respect 
to L. That is, f(x) follows CAA(0, 1). For this case, it can be shown that E[| det(Hi)|] = ^ and, thus, the theoretical 
limit in Theorem [5] leads 4 — 41og(% L ) (approximately 3.1). The detailed proof of E[| det(Hi)|] = ^ is in Appendix 
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II. As shown in the figure, i? m i mo — Rm quickly converges to the theoretical limit as L increases. Considering that 
the sum capacity gap is approximately given by 4.7 bits/sec/Hz when L = 2 (Theorem [3] and Example [3]), the sum 
capacity gap can be tightened as L increases. 

D. Approximate Ergodic Capacity for Fading Interference Channel 

We notice that a similar analysis used in Theorems [2] and [3] is applicable to show an approximate ergodic 
capacity for fading K-user interference channel. The achievability follows from ergodic interference alignment in 
|[3l . Assuming that all sources employ uniform power allocation across time, we show that ergodic interference 
alignment characterizes an approximate ergodic per-user capacity, i.e., ergodic sum capacity divided by K, for 
a broad class of channel distributions. The detailed statement is given in Theorem [6] in Appendix III. For i.i.d. 
Rayleigh fading, for instance, our analysis characterizes the ergodic per-user capacity to within ^ log 6 bits/sec/Hz 
(approximately 1.3 bits/sec/Hz). 



IV. Ergodic Interference Neutralization 

For the achievability, we propose ergodic interference neutralization using an even number of relays. Let M := 
L|j.Then we can choose 2M relays among the total number L of relays and apply the proposed ergodic interference 
neutralization by using these 2M relays. For simplicity, we assume L is even in the rest of this section. That is, 
L = 2M. 



A. High-Level View 

Before the detailed description and analysis, we begin by providing a high-level view of the proposed ergodic 
interference neutralization. Consider length-n sequences of matrices {H[t]}" =1 and {G[i]}™ =1 , drawn i.i.d. ac- 
cording to a certain probability density function. We partition these sequences judiciously into pairs of matrices 
(H[ii], Gfo]) such that Gfo] and F(H{ti]) are almost equal, where F(-) is a cleverly chosen mapping to be 
discussed below. The main argument is that by considering a longer and longer sequence of matrices, we can make 
these two matrices arbitrarily close. The formal and technical details of this argument can be found in Sections 
HV-Bl and IIV-CI For notational convenience, we introduce the notation G[t 2 ] ^ F(H[tj]) for the two matrices that 
are almost equal. 

As pointed out in [6], a simple amplify-and-forward scheme with an appropriate delay r £ Z + can neutralize 
interference by letting G[t + r]H[t] approximately a diagonal matrix with non-zero diagonal elements. To satisfy 
this condition, we first partition L relays into M = 4 pairs and neutralize interference separately by each pair of 
relays. Figure [6] illustrates the main idea of the proposed scheme. For A = {aj j} € C 2x2 , define 



Fo(A) 



0-2,2 «1,2 



(15) 
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The relays then amplify and forward with delay r such that G m [t + r] ~ i<2(H m [i]) for all m G {1, • • • , M}. For 
relaying, relays 2m — 1 and 2m amplify and forward with the amplification factors 7 ^t^H* Tq | an ^ ~7 ra^g n? ] | ■ 

respectively. Here 7 = \J 1 + 2 p * s nee ded to satisfy the average power constraint P. Then the effective channel 
matrix of the mth pair is given by 



. >t(H m [t])* 
I det(H m [t])| 



1 

-1 



1 
-1 



H m [t] 



7l det(H r 



1 
-1 



(16) 



As a consequence, the effective channel gain from each source to its destination is approximately given by 
T 2 (£m=i I det(H m [i])|) 2 , as can be seen in the rate expression in Theorem Q] One can easily show that the 
additional noise power at destination i due to this amplify-and-forward relaying is given as a\ ¥ p as shown in the 
rate expression in Theorem [T] Lastly, since the probability density functions of the paired channel states are the 
same, i.e., 

/ H[ t]([Hf,--- ,Hl I ] T ) = f G[t] ([F 2 (U 1 ),--- ,F 2 (H M )]), (17) 

almost all channel instances can be utilized for this ergodic pairing as the block length n increases. Hence, the 
ergodic rate in Theorem Q] is achievable in the limit of large n. 

There are two crucial facts to be observed: 1) the intended signal power received at each destination is non-zero 
while the interference power decreases arbitrarily close to zero at any finite power P; 2) the intended signal power 
received at each destination increases quadratically with increasing L. These facts make approximate capacity 
characterization possible for a broad class of channel distributions. 

Although finding a pair of channel instances having exact prescribed values is impossible, such a pairing can be 
done approximately by partitioning the channel space of each hop and then pairing the partitioned channel spaces 
between the first and second hops. In the following subsection, we first explain channel space partition and pairing 
and then explain the detailed scheme. 

B. Block-Wise Ergodic Interference Neutralization 

1) Partitioning and pairing of channel space: We partition the channel space of each hop, i.e, C 2Mx2 space for 
the first hop and C 2x2A/ space for the second hop. First, consider the channel space of the first hop C 2A/x2 . For 
N G Z + and A > 0, define 

Ql :={A G A(Z 2Mx2 + jZ 2Mx2 )||re(a ij -)| < AN, \im(a id )\ < AN 

for all i G {1, • • • ,2M} and j G {1,2}}, (18) 

where A = {ojj}. Here, N and A are related to the number of quantization points and the quantization interval. 
For a quantized channel matrix Q G Q\, define 

-4i(Q) :=|a g C 2Mx2 | - y < re(oy) -re(q M ) < ^ and - y < im(aij) - im^y) < y 

for all i G {1, • • • , 2M} and j G {1, 2} 1, (19) 



where A = {ctij} and Q = {qij}- Figure [7J illustrates the channel space partitioning with respect to hij G C. We 
can define Q2 and .42 (Q) f° r the second hop as the same manner in (ED and C[9]) by substituting A G A(Z 2x2M + 
j7? x2M ) and A G C 2x2M , respectively. We will only use the first-hop channel instances in Uq 6 q 1 ^4i(Q) and the 
second-hop channel instances in UQ g g 2 ^l2(Q) for transmission. 

Now consider the channel space pairing between A\(Q) and ^(Q). For A G C 2Mx2 , define 

F{A) := [F 2 (A 1 ),F 2 (A 2 ),--- ,F 2 (A M )], (20) 

where A = [Af , Aj, • • • , A^ f ] T and the definition of F 2 (-) is given by (fT5T ). For H[i] G A\(Q), the relays will 
amplify and forward with delay r G 1* + satisfying G[t + r] G ^4 2 (F(Q)). Hence the channel subspace ^4i(Q) of 
the first hop is paired with the channel subspace „4 2 (F(Q)) of the second hop. The detailed transmission scheme 
is given in the following subsection. 
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AN 



QeSi 
-A(Q) 



re(hi 



Fig. 7. Channel space partitioning with respect to the channel coefficient hi : j £ 



2) Transmission scheme: We first divide a length-n block into B sub-blocks having length rig = each. At 
the first sub-block, the sources transmit their first messages to the relays (the relays do not transmit). At the 6th 
sub-block, b G {2, • • • , B — 1}, the sources transmit their 6th messages to the relays and the relays amplify and 
forward the received signals of the (6— l)th sub-block to the destinations. At the last sub-block, the relays amplify 
and forward the received signals of the (B — l)th sub-block to the destinations (the sources do not transmit). Hence, 
the number of effective sub-blocks is equal to B — 1. Since we can set both ng and B as large as desired as n 
increases, the fractional rate loss becomes negligible as n increases. For simplicity, we describe the proposed 
scheme based on the first message transmission and omit the sub-block index. 

For Q € Qi, define 71(Q) := {t G {1, • • • , n,e}|H[i] G ^4i(Q)}, which is the set of time indices of the first hop 
whose channel instances belong to „4i(Q). Similarly, for Q G Q 2 , 7i(Q) := {t 6 {ns + 1, • • • , 2ns}|G[t] G ^(Q)}, 
which is the set of time indices of the second hop whose channel instances belong to ^(Q). The encoding, relaying, 
and decoding are as follows. 

• (Encoding) The sources transmit their messages using Gaussian codebook with length n b and average power 
P. 

• (Relaying) For all Q G Qi, the relays amplify and forward their received signals that were received during 
71(Q) using the time indices in ^(^(Q))- Specifically, for t\ € 71(Q), the transmit signal vector of the 
relays is given by x fl [t 2 ] = ry^ti], where i 2 G ^(^(Q))- Here 

,jgt(Q0! a 



7 |det(Q 1 )|- 








dct(Q 2 )* 
7 |det(Q 2 )|" 







dct(Q M )* a 
T|det(Q M )r V 



(21) 



, Q^] T and denotes the 2 x 2 dimensional 



7 = y j^p, and A = [[1, 0} 1 [0, -I} 1 Y , where Q = [Q^ 
all-zero matrix. 

(Decoding) The destinations decode their messages based on their received signals during t G {n^+1, • • • , 2n#}. 



C. Achievable Rate Region 

In this subsection, we prove Theorem [TJ We first introduce the following two lemmas. 

Lemma 1: For any Q G Qi, 

P[H[i]GA(Q)] = P[G[i]G^ 2 (F(Q))]. 



(22) 
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Proof: Let fu[t]{') an d /g[*](") denote the probability density functions of H[t] and G[t], respectively. Then 

P[H[t]€Ai(Q)]= / /H[*](A)dA 
•/Ae.A 1 fQl 



/ II /(aij)dA 

./Ae.4i(Q) ie {i i ... 2 Af}je{l,2} 



= / / G[t] (F(A))dA 

= / /G N (A')dA' 
iA'eA(F(Q)) 

= P [G[t] G -4(F(Q))] , (23) 

where A = {dij}. Here (a) follows from the definition of F(A) and (6) follows by a change of variable A' = F(A) 
whose Jacobian is one and „4 2 (F(Q)) = {F(A)\A G Ai(Q)}. Therefore Lemma Q] holds. ■ 



Lemma 2: The probability that 

card(Ti(Qi)) 



P[U[t] G -4_i(Qi)] 

n B 



and 

card(r 2 (Q 2 )) 



P[G[t] g .4 2 (Q 2 ) 



< S (24) 



< 5 (25) 



n B 

for all Qi G Qi and Q 2 G Q 2 is greater than 1 - (card(Qi) + card(Q 2 ))/(2n B 5 2 ). 

Proof: We refer to Lemma 2.12 in [30 ] for the proof. ■ 

Suppose that the sources transmit at time t\ G 71 (Q) and the relays amplify and forward their received signals at 
time t 2 G T 2 (F(Q)), where Q G Q\. For this case, from (Q]) and (O, the received signal vector of the destinations 
is given by 

y[fe] = G[t2]rH[ti]x[*i] + G[t 2 ]rz fl [t 1 ]+z[t 2 ], (26) 

where we use x R [t 2 ] = Ty R [h]. Denote H[ti] = H = [Hf,---H^] T and G[t 2 ] = F(H) + A, where A = 
[Ai, • • • , Am] is the quantization error matrix with respect to F(H). From (1261) . 

y[*a] = ( (l J2 I det ( H ™)l) A + ArH ) X N + ( F ( H ) + A ) rz i?N + Z N, (27) 

where we use F(H)TH = (7 2~2m=l I det(H m )|)A. Thus, the received signal-to-interference-and-noise ratio (SINR) 
of destination i is given by 

■l) i - 1 (7E^=i|det(H m )|) + [ArH] ; 



P 

SINR, 



(28) 



1 + 7 =1 (I [H m ]3-i,3-i + [A m ] M | 2 + I [H m ] ii3 _i + [A TO ] i)3 _i| 2 ) + P| [ArH]i )3 _i| 2 ' 
Define Ri(Q) = minAe,A 1 (Q) log(l + SINRj). Then an achievable rate of destination i is lower bounded by 

Ri > — J2 i?i(Q)min{card(Ti(Q)),card(r 2 (F(Q)))}. (29) 
QeQi 

From Lemmas [J and |2l 

card(Ti(Q)) > n B (P[H[t] G A(Q)] - 5) (30) 

and 

card(r 2 (F(Q))) > n B (P[G[t] G A 2 (F(Q))} - 5) 

= n B {P[H[t]eA l {Q)}-8) (31) 
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for all Q S Qi with probability greater than 1 — — ■ wnere we use card(Qi) = card(Q2) = (2N + 1) 8M . 

Then 



Ri > Yl ^(Q)(P[HM g Ai(Q)} - 6) 



QeQi 

> Y ^(Q) P [ H W G ^(Q)] - S2 ( 2N + l) 8M maxi2i(Q) 
QeQi QeSl 



is achievable with probability greater than 1 

— 1/3 

5 = n B , the following condition can be satisfied: 

A = n 



(2N+1) 8 



n B 8 2 



By setting A 



-1/(3- 2 5 M) 
B 



n 



l/(3-2 4 M) 
B 



(32) 



, and 



AiV = n 



-l/(3-2 5 M) 
l/(3-2 5 M) 



OO, 



62(2N + lf M max RAQ) < 2 • 3 8M JV 8M <5 max JUQ) 
QeQi QeQi 

(a) 



(2iV + l) 8M < s 8M N 8M 



< 2 • 3 8M iy 8M <51og(l + 2 4 MA 2 iV 2 P) 
2 • 3 8M n B 1/6 log(l + 2 4 M n ; /(3 ' 24M) P) -> 0, 
3 8M n -i/6 ^ Q 



(33) 



ns<5 2 ne<5 2 

as ris increases, where (a) follows since \hij\ 2 < 2A 2 (N + |) 2 < 2 3 A 2 iV 2 for the channel instances using the 
transmission (see Fig. [7]). 
Hence, 



Ri 



( 



log 



1 + 



Pi 2 E^=i|det(H r 



\ 



V 



1 + 7 2 E^=l(l^2(m-l)+3-i,3-i| 2 + l^2(m-l)+i,3-i| 2 ) 

is achievable with probability approaching one for i G {1, 2}, where we use the fact that 

Pi 2 feLi|det(H m )|) 2 
lim SINRj = Wl ^ '- . 

1 + 7 2 E m =l(l^2(m-l)+3- 4 ,3-i| 2 + \ h 2( m ~l)+i,3-i\ 2 ) 

In conclusion, Theorem Q] holds. 



(34) 



(35) 



V. Approximate Capacity Characterization 

In this section, we prove Theorems [2] to [5J the approximate ergodic sum capacity characterization results. We 
will deal with the difference between i? m i mo and R m , which are given by © and (O respectively. Throughout 
this section, we assume a class of channel distributions such that f(x) is only a function of \x\. That is, for given 
amplitudes of the channel coefficients, their phases are i.i.d. uniformly distributed over [0, 2tt). For instance, this 
class of channel distributions includes i.i.d. uniform phase fading and i.i.d. Rayleigh fading as special cases. We 
omit the time index t in this section for notational convenience. 



A. Approximate Capacity for L = 2 

We first consider the case where L = 2. In order to deal with i.i.d. random phase in the rate expression in 
Theorem [TJ we introduce the following lemma showing the exact solution of [log (1 — xcosc/))] for \x\ < 1 
when (j) is uniformly distributed over [0, 2?r). 

Lemma 3: Let be a random variable uniformly distributed over [0, 27r). For \x\ < 1, 

E [log (1 - x cos <p)} = log (l + y/l - x 2 ) - 1. (36) 



Proof: We refer to the equation (4.224 12) in EH. 
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1) Proof of Theorem^ From ([8]), 



R m = 2 Efi 



log 1 + 



2P 2 (1 -cos0) 



2 log 1 + 



2P 2 



2 log 1 + 



1 + 4P 

2P 2 
1 + 4P 



1 + 4P 



2E fi 



log 1 



2P 2 



1 + 4P + 2P 2 



cos 9 



+ 21og 1 + Wl 



2P 2 



1 + 4P + 2P 2 



(37) 



where 9 = 9\ t i + 6> 2 ,2 — 6*1,2 — #2,1- Here, (a) follows since | det(H)| 2 = 2(1 — cos 6) and <r^ Fi = 1 ^ p , (b) follows 
since 6 mod[27r] is uniformly distributed over [0, 2ir) and from Lemma [3] with | \ + 4P+2P 2 1 — ^ Similarly, from 



Pn 



E e [log((l + 2P) 2 - 2P 2 (1 + cos 6))] 



log(l + 4P + 2P 2 ) + E e 



log 1 



2P 2 



1 + 4P + 2P 2 



cos 9 



log(l + 4P + 2P 2 ) + log 1 + 1/1 



2P 2 



1 + 4P + 2P 2 



(38) 



Then, from <|37> and (f38]), 



log 



(a) 

< log 
(*>) 

< 4, 

2P : 



(1+4P) 2 
1 + 4P + 2P 2 

(1 + 4P) 2 
1 + 4P + 2P 2 



log 1 + 




2P 2 



1 + 4P + 2P 2 



+ 1 



+ 1 



(39) 



where (a) follows since | i + 4P+2P 2 1 — ^ ^ or an ^ ^ > and (b) follows since 



log 



(1 + 4P) 2 
1 + 4P + 2P 2 



< log 

= 2 log 

< 3, 



(1 + 4P) 5 



1 + 2\/2P + 2P 2 
1 + 4P \ 

1 + y/2Pj 



(40) 



where we use the fact that log (^ j^j^p ) * s an increasing function of P > and linip^oo log ( ^^p ^j = f • m 
conclusion, Theorem [2] holds. 

2) Proof of Theorem\3} Since f(x) is only a function of \x\, hij can be represented as Ojj exp^fljj), where 
fljj > and G [0, 27r) are independent of each other. Moreover 9^j is uniformly distributed over [0, 27r). 
To simplify the notation, we denote a = {01.1,01,2,02,1,02.2}, A = a\ l a\ 2 + af 2 a 2i> = a i 1 + a l 1 + 
P 2 = af 2 + a | 2 + 2 > C = 201,101,202,102,2, and 5 = of x + of 2 + a 2 , x + a| 2 . 
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From ®, 



(a) 



(b) 



E e. 

ie{l,2} 

E E - 

ie(l,2} 

E e. 

'•6(1.2} 



log 1 



P 2 (A-Gcos9) 



l + PBi 



log 1 + 



log 1 + 



^ E E * 

i6(l,2} 

= E ^8 

16(1,2} 



log ^1 + 
log ( 1 + 



P 2 A 
l + PBi 

P 2 A 
l + PBi 

P 2 A 
l + PB, 

PA\ 
~B~) 



+ E E * 

16(1,2} 

+ E E * 

16(1,2} 



log 1 



P 2 Gcosi 



1 + PBi + P 2 a) 



+ 2E 8 



log U + 



log I i + yi- 

VA 2 -G 2 



P 2 G 



A 



1 + PB l + P 2 A 
- 2 



+ E Ea 

i6(l,2} 



log 



+ 2E a 



log 1 + 



VA 2 - G 2 



> E E * 

ie{l,2} 



log 1 + 



4 
PM\ 



+ 2E E 



log 1 + 



VA 2 - G 2 



+ E Ea 

i6{l,2} 



log 



Bi + P(A + B 2 ) + P 2 ABi 



(41) 



where 6 = + #2,2 — #1,2 — #2,1 ■ Here (a) follows from the facts that a and {#1,1, #1,2, #2,1, ^2,2} are independent 
of each other and | det(H)| 2 = A — Gcos 0, (6) follows since 6 mod [27r] is uniformly distributed over [0, 2it) and 
from Lemma [3] with l+p P g ( ^p2 A < 1. (c) follows since — 

pbUp^a ^ S for 

any P > 0, and (d) follows since 

From ®, 



> log 77 for ci, c 2 , c 3 > and c 2 < c 3 . 



R ( - } F 



logdet^ + PHH"!) 



= E a [Ee [log (1 + PS + P 2 (A - Gcosd))]] 



= E a [log(l + PS + P 2 A)] + E a 
( = } E a [log(l + PS + P 2 A)} + E a 
< E a [log(l + P5 + P 2 ^)] , 



log 1 



P 2 G 



1 + PS + P 2 A 



cos 6)) 



log 1 + 4 1 



P 2 G V 

1 + ps + p 2 ^ y 



(42) 



where (a) follows from the fact that a and {#1,1, #1,2, #2, 1, 6*2.2} are independent of each other, (6) follows since 
9 mod[27r] is uniformly distributed over [0, 27r) and from Lemma [3] and 1 i+ps+p2A I < 1, and (c) follows again 



since 



P 2 G 



Let 



1+P5+P 2 A 



< 1. 



A = log(l + PS + P 2 ,4)- E lo s( 1 + ^)- 



(43) 



i6(l,2} 
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Then 



A loo [ ^1 ) + ] og 



< log 



A 
A 



1 + PS + P 2 A 



^ + P(B l + B 2 ) + P 2 A 



where the inequality follows since P>iE>2 > A and B\ + B 2 > S. Therefore, from (|4TT > to (l44l . 

£? 2 



-Rmimo — -Rin < E a [A] — ^ E, 

i6{l,2} 



log 



2E, 



log 
log 



A + B 2 

A{A + Bl){A + Bl) 
B^A + VA^-G 2 ) 2 

f Va(a + b 2 ) \ 
[b^a + va^-g 2 ) 



-2E 8 

+ 2 
+ 2. 



log 1 + 



VA 2 - G 2 



+ 2 



In conclusion, Theorem [3] holds. 



(44) 



(45) 



B. Approximate Capacity as L — > oo 

In this subsection, we characterize an approximate ergodic sum capacity in the limit of large number of relays by 
deriving lirn/^ooji^mimo — i?i n }. For .fT-user 2-hop networks with L relays, it was shown in O that interference 
can be completely neutralized if K > N(N — 1) + 1, which indicates that for 2 x L x 2 networks interference 
neutralization can be achieved without channel pairing if L > 3. However, maximizing the achievable sum rate 
exploiting interference neutralization without channel pairing presented in |[32l is non-convex and, as a result, 
it is unclear how to determine the sum rate gap from the cut-set upper bound. By contrast, we now show that 
our achievable rate expression from Theorem Q] permits to derive a finite-gap result. The rate expression R m in 
dD contains the sum of i.i.d. random variables, i.e., Y^m=i |det(H m )|, which approaches a deterministic value 
M E[| det(Hi)|] almost surely as M — > oo by the law of large numbers. The following lemma provides a rigorous 
lower bound in order to deal with i?j n that holds for any M. 

Lemma 4: Consider a sequence of i.i.d. nonnegative random variables {Xi,i G ^+}- Let S m = Y^ILi^i- If 
E[Xi 2 ] < oo, then for any e G (0, E[Xl]) and any c > 0, 

E [log(l + cS m 2 )] > log (1 + cm 2 (E[X!]) 2 ) - 5 m (c, E[Xi], E[Xi 2 ]), (46) 

where 

Mc, E[*i], E[X X 2 ]) = log (1 + C m 2 (E[Xx] - e ) 2 ) 

is a positive sequence of m, which converges to zero as e — > 0. 
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Proof: We have 



E [log(l + cS m 2 )] 

= E [log (l + cS m 2 ) (l{|S m /m-E[Xx]|<6} + 1 {|5 m /m-E[X 1 ]|>e})] 

> log (1 + C m 2 (E[X!]) 2 - 2cm 2 e E[X X ] + cm 2 e 2 )) E[l { | 5m/m _ E[Xl] | <e} ] 

> log (1 + cm^EfXx]) 2 - 2cm 2 e E[X X ] + cm 2 e 2 )) f 1 - Y^li 

\ me* 

(c) 

> log (1 + cm 2 (E[Xi]) 2 - 2cm 2 e E[XJ + cm 2 e 2 )) 
, log(l + cm 2 (E[X 1 ]) 2 + cm 2 e 2 )) 



me 4 



log (1 + cm 2 (E[X!]) 2 ) - <5 m (c, E[Xi], E[*i 2 ]), 



(48) 



where (a) follows since 5m > m E[X\] —me under the condition \S m /m— E[Xi]| < e, (6) follows from Chebyshev's 
inequality, and (c) follows since Var(Xi) < E[Xi 2 ]. In conclusion, Lemma [4] holds. ■ 
By setting S m arbitrarily small as m increases, Lemma |4] provides that 

E [log(l + cS m 2 )] > log (1 + cm 2 (E[X!]) 2 ) (49) 

in the limit of large m. Note that this bound is asymptotically tight since E [log(l + cS m 2 )] < log(l + cE[5 m 2 ]) 
from Jensen's inequality and log(l + cE[5 m 2 ]) is approximately given as log(l + cm 2 E[Xi] 2 ) as m increases. 



1) Proof of Theorem® Recall M = That is, L < 2M + 1. From ©, 

Pmimo < 21og(l + P(2Af +1)), 



(50) 



where the inequality follows from Jensen's inequality and the fact that logdet(-) is a concave function |f33l . Here 
we assume L = 2M + 1 to obtain an upper bound. 
From ®, 



(a) 



2E 



{0i 



log 



1 + 



V 



P 2 (E^i V2^2^T, 
1 + P(2M + 2) 



(*>) / 
> 21og 1 + 



Wp2 M 2 



1 + P(2M + 2) 



25 M 



P 1 



1 + P(2M + 2)'tt 



(51) 



where 9 m = 02m-i,i + 02m,2 — #2m-i,2 _ #2m,i- Here, (a) follows since | det(H m )| = \/2 — 2 cos 6 m and (6) follows 
since 6 m mod [27r] is uniformly distributed over [0, 27r) and from Lemma [4] with the facts that E [y/2 — 2 cos 9\] = 
and E [2 - 2 cos 6>i] = 2. Then, from (gOj and (gTj), 



Pn 



P in < 2 loe 



+ 2,5 M 



1 + P(4M + 3) + P 2 (2M + 1)(2M + 2) 
1 + P(2M + 2) + i|P 2 M 2 
P 2 4 



Hence, liniM^oo{Pn 



l + P(2M + 2)'vr : 
Pin} < 4 log 7T — 4 + e\, where 

ei= i i g 1 oo^ M ( l + P(2M + 2) '^ 

.2 

' 1 >o, 



(52) 



-2 log 1 - e 



7T 7T 



2 16 

which can be arbitrarily small as e decreases. In conclusion, Theorem @] holds. 



(53) 
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2) Proof of Theorem\5} From ([8]), 

R m > 2 E 



log L + P 2 ^|det(H m )|j 

/ / 2M \ 

log 1 + P [J2 \hm,i\ 2 + 2 



\m=l 



( >21og(- + jP M (E[|d6t(Hl ' 



1 + P(2M + 2) J 
-26 M (P 2 ,E[\ detCHOH.EjldetCHOI 
where (a) follows from Lemma [4] and Jensen's inequality. Hence, from (TSUI ) and (f54T >. 

'1 + P(4M + 3) + P 2 (2M + 1)(2M + 2) 



(54) 



R 



mimo x '-in 



Pi, < 2 log 



1 + P 2 M 2 (E[ 


detlHOI]) 2 


(p 2 ,E[|det(H 1 )|],E 


JdetCHOI 2 ]) 



and 



where 



hill {Pmimo 

M->oo 



Pin} <4-41og(E[|det(Hi)|]) + e 2 , 



(55) 



(56) 



e 2 = lim 25 M (p 2 ,E[|det(H 1 )|],E[|det(H 1 )| 2 l) 

M— >oo \ L J / 

21 °H (E[|det( Hl )|]) 2 ; >U ' 

which can be arbitrarily small as e decreases. In conclusion, Theorem [5] holds. 



(57) 



VI. Conclusion 

In this paper, we studied a fading 2-user 2-hop network with L relays where channel coefficients vary over 
time. In spite of recent achievements in this area, the best known capacity characterization is to within o(log SNR) 
bits/sec/Hz from the ergodic sum capacity, which can be arbitrarily large as SNR increases. For a broad class 
of channel distributions, we tightened this gap to within a constant number of bits/sec/Hz, independent of SNR. 
The achievability follows from ergodic interference neutralization in which the relays are partitioned into several 
pairs and interference is neutralized separately by each pair of relays. The proposed scheme makes interference 
neutralized in the finite SNR regime and, at the same time, the intended signal power increased quadratically with 
L, leading that the optimal 21og(LSNR) rate scaling is achievable, which cannot be captured by the previous DoF 
work. 

Appendix I 

Quantization for i.i.d. Uniform Phase Fading 

For i.i.d. uniform phase fading, hij[t] and gj^t] respectively are represented as exjp(j6ij[t]) and exp(jtpj t i[t]) 
for i G {!,■■■ ,L} and j G {1,2}. Hence we can quantize the channel space of each hop based on angles. 
Specifically, the channel space of the first hop can be partitioned as follows. For N G Z+, first define Q\ := 
{ exp(jO), exp(j^), exp(j^), • • • , exjp(j ^ N ~^ 27T )} 2 x2 . Let U denote the set of all x G C satisfying \x\ = 1. 
For a quantized channel matrix Q G Q\, define ^li(Q) := {A G U 2Mx2 \ — j; < ^a-ij - ^q%,j < 77 for all i G 
{1, • • • , 2M} and j G {1,2}}, where A = {aj j}, Q = {qij}, and Zx denotes the angle of x G U, i.e., x = 
exp(jZx). Figure [8] illustrates the channel space partitioning with respect to hij G U. In a similar manner, we can 
define Q2 and .42 (Q) f° r the second hop. Then we can show that there exists an increasing sequence of N, which 
is a function of rig, such that (f7]) is achievable as increases using similar steps in the proof of Theorem [T] 
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Fig. 8. Channel space partitioning with respect to the channel coefficient hij £ U for i.i.d. uniform phase fading. 



Appendix II 

Closed Form of E[| det(Hi)|] for i.i.d. Rayleigh Fading 

Let A be a 2 x 2 matrix whose entries are i.i.d. circularly symmetric complex Gaussian random variables with 
mean zero and unit variance and W := 2AA"I\ Let Ai and A2, Ai > A2, be the eigenvalues of W. Then the joint 
probability density function of Ai and A2 is given by 041 Equation (3.11)] 



Thus, 



/(Ai, A 2 ) = ^ exp + A 2 )^ (A x - A 2 ) 2 l {Al > A2 > 0} (A!, A 2 ). 
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(58) 



(59) 



where (a) follows by a change of variable u = (Ai + A 2 )/2 and v = (Ai — A 2 )/2. In conclusion, E[| det(H)|] = ^ 
for i.i.d. Rayleigh fading. 



Appendix III 

Approximate Ergodic Capacity for Fading Interference Channel 

A similar analysis used in Theorems [2] and [3] is applicable for fading interference channel. Specifically, consider 
the K-user interference channel in which the input-output relation is given by 



y[t]=H[t]x[t]+z[t] 



(60) 
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and the elements of H[i] = are i.i.d. drawn from a continuous distribution f(x), x G C, and vary 

independently over time. The achievability follows from ergodic interference alignment in [3] showing that 

1 



Ri 



E[log(l + 2\httfP)] 



(61) 



is achievable for all i G {1, • • • , if} (3] Theorem 2]. Theorem [6] characterizes an approximate ergodic per-user 
capacity, i.e., ergodic sum capacity divided by K, assuming that all sources employ uniform power allocation across 
time. For this case, the sum of any pair of achievable rates is upper bounded by 



Ri + Rj < E 



log 1 + 



mm 



(62) 



for all i,j G {1, • • ■ ,K}, i ^ j |3] Equation (99)]. From the lower bound (loTT ) and the upper bound (l62l . we 
characterize an approximate ergodic per-user capacity in the following theorem. 

Theorem 6: Consider the fading K-user interference channel. Let R 1SL := Yl%=i \ Epog(l + 2\hi^\ 2 P)} and C sum 
denote the sum capacity assuming that all sources employ uniform power allocation across time. Then 
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for any P > 0. 
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Then, from (|6Tb and (l62l and the fact that channel coefficients are i.i.d., 
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The first term of (1661) is upper bounded as 
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Similarly, the second term of d66l ) is upper bounded as 
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where (a) follows since amax{l,6/a} = 6max{l,a/6} for all a, b > 0. Therefore, 
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Finally, combining (1651 ) and (1691 ) shows the gap in (1631 ). which completes the proof. 
For i.i.d. Rayleigh fading channel, |/iij| 2 has the exponential distribution and 
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and the gap in Theorem [6] is given as ilog6 bits/sec/Hz (approximately 1.3 bits/sec/Hz). 



(68) 



(69) 



(70) 



References 

[1] R. H. Etkin, D. N. C. Tse, and H. Wang, "Gaussian interference channel capacity to within one bit," IEEE Trans. Inf. Theory, vol. 54, 
pp. 5534-5562, Dec. 2008. 

[2] V. R. Cadambe and S. A. Jafar, "Interference alignment and degrees of freedom of the A"-user interference channel," IEEE Trans. Inf. 

Theory, vol. 54, pp. 3425-3441, Aug. 2008. 
[3] B. Nazer, M. Gastpar, S. A. Jafar, and S. Vishwanath, "Ergodic interference alignment," IEEE Trans. Inf. Theory, vol. 58, pp. 6355-6371, 

Oct. 2012. 

[4] S. Mohajer, S. Diggavi, C. Fragouli, and D. N. C. Tse, "Approximate capacity of a class of Gaussian interference-relay networks," 

IEEE Trans. Inf. Theory, vol. 57, pp. 2837-2864, May 2011. 
[5] T. Gou, S. A. Jafar, C. Wang, S.-W. Jeon, and S.-Y. Chung, "Aligned interference neutralization and the degrees of freedom of the 

2x2x2 interference channel," IEEE Trans. Inf. Theory, vol. 58, pp. 4381^1395, Jul. 2012. 



22 



[6] S.-W. Jeon, S.-Y. Chung, and S. A. Jafar, "Degrees of freedom region of a class of multisource Gaussian relay networks," IEEE Trans. 

Inf. Theory, vol. 57, pp. 3032-3044, May 2011. 
[7] C. Wang, T. Gou, and S. A. Jafar, "Multiple unicast capacity of 2-source 2-sink networks," in arXiv:cs. IT/1 104.0954, Apr. 2011. 
[8] I. Shomorony and A. S. Avestimehr, "Two-unicast wireless networks: Characterizing the degrees-of-freedom," in arXiv:cs. IT/1 102.2498, 

Feb. 2011. 

[9] B. Rankov and A. Wittneben, "Spectral efficient protocols for half-duplex fading relay channels," IEEE J. Sel. Areas Commun., vol. 25, 
pp. 379-389, Feb. 2007. 

[10] V. R. Cadambe and S. A. Jafar, "Interference alignment and the degrees of freedom of wireless X networks," IEEE Trans. Inf. Theory, 

vol. 55, pp. 3893-3908, Sep. 2009. 
[11] B. Nazer, M. Gastpar, S. A. Jafar, and S. Vishwanath, "Ergodic interference alignment," in Proc. IEEE Int. Symp. Information Theory 

(ISIT), Seoul, South Korea, Jun./Jul. 2009. 
[12] M. A. Maddah-Ali, A. S. Motahari, and A. K. Khandani, "Communication over MIMO X channels: Interference alignment, 

decomposition, and performance analysis," IEEE Trans. Inf. Theory, vol. 54, pp. 3457-3470, Aug. 2008. 
[13] T. Gou and S. A. Jafar, "Degrees of freedom of the K user M x N MIMO interference channel," IEEE Trans. Inf. Theory, vol. 56, 

pp. 6040-6057, Dec. 2010. 

[14] C. Suh, M. Ho, and D. N. C. Tse, "Downlink interference alignment," IEEE Trans. Commun., vol. 59, pp. 2616-2626, Sep. 2011. 
[15] C. Suh and D. N. C. Tse, "Interference alignment for cellular networks," in Proc. 46th Anna. Allerton Conf. Communication, Control, 

and Computing, Monticello, IL, Sep. 2008. 
[16] V. R. Cadambe and S. A. Jafar, "Degrees of freedom of wireless networks with relays, feedback, cooperation, and full duplex operation," 

IEEE Trans. Inf. Theory, vol. 55, pp. 2334-2344, May 2009. 
[17] V. S. Annapureddy, A. El Gamal, and V. V. Veeravalli, "Degrees of freedom of interference channels with CoMP transmission and 

reception," IEEE Trans. Inf. Theory, vol. 58, pp. 5740-5760, Sep. 2012. 
[18] L. Ke, A. Ramamoorthy, Z. Wang, and H. Yin, "Degrees of freedom region for an interference network with general message demands," 

IEEE Trans. Inf. Theory, vol. 58, pp. 3787-3797, Jun. 2012. 
[19] A. S. Motahari, S. O. Gharan, and A. K. Khandani, "Real interference alignment with real numbers," in arXiv:cs. IT/0908. 1208, 2009. 
[20] A. S. Motahari, S. O. Gharan, M. A. Maddah-Ali, and A. K. Khandani, "Real interference alignment: Exploiting the potential of single 

antenna systems," in arXiv.cs.1T/0908.2282, 2009. 
[21] G. Bresler, A. Parekh, and D. N. C. Tse, "The approximate capacity of the many-to-one and one-to-many Gaussian interference 

channels," IEEE Trans. Inf. Theory, vol. 56, pp. 4566^1592, Sep. 2010. 
[22] W. Nam, S.-Y. Chung, and Y. H. Lee, "Capacity of the Gaussian two-way relay channel to within | bit," IEEE Trans. Inf. Theory, 

vol. 56, pp. 5488-5494, Nov. 2010. 
[23] A. S. Avestimehr, S. N. Diggavi, and D. N. C. Tse, "Wireless network information flow: A deterministic approach," IEEE Trans. Inf. 

Theory, vol. 57, pp. 1872-1905, Apr. 2011. 
[24] S. H. Lim, Y.-H. Kim, A. El Gamal, and S.-Y. Chung, "Noisy network coding," IEEE Trans. Inf. Theory, vol. 57, pp. 3132 - 3152, 

May 2011. 

[25] U. Niesen and M. A. Maddah-Ali, "Interference alignment: From degrees-of-freedom to constant-gap capacity approximations," in 
arXiv.cs.IT/l 112.4879, 2011. 

[26] O. Ordentlich, U. Erez, and B. Nazer, "The approximate sum capacity of the symmetric Gaussian fc-user interference channel," in 
arXiv.cs.IT/1206.0197, 2012. 

[27] S.-W. Jeon and S.-Y. Chung, "Capacity of a class of multi-source relay networks," in Information Theory and Applications Workshop, 

University of California San Diego, La Jolla , CA, Feb. 2009. 
[28] U. Niesen, B. Nazer, and P. Whiting, "Computation alignment: Capacity approximation without noise accumulation," in 

arXiv:cs.lT/l 108.6312, 2011. 

[29] I. E. Telatar, "Capacity of multi-antenna Gaussian channels," European Trans, on Telecommun., vol. 10, pp. 585-595, Nov. 1999. 
[30] I. Csiszar and J. Korner, Information Theory: Coding Theorems for Discrete Memoryless Systems. New York: Academic Press, 1981. 
[31] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products. Academic Press, 2007. 

[32] C. Esli, S. Berger, and A. Wittneben, "Optimizing zero-forcing based gain allocation for wireless multiuser networks," in Proc. IEEE 

International Conference on Communications (ICC), Beijing, China, Jun. 2007. 
[33] S. Boyd and L. Vandenberghe, Convex Optimization. New York: Cambridge Univ. Press, 2004. 

[34] A. Edelman, "Eigenvalues and condition numbers of random matrices," Ph.D. dissertation, Massachusetts Institute of Technology (MIT), 
1989. 



